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£Sj ■ We point out that ferromagnetic superconductors with triplet pairing and strong spin-orbit cou- 

' pling are even in the simplest case at least two-band superconductors. The Gor'kov type formalism 

for such superconductors is developed and the Ginzburg-Landau equations are derived. The depen- 
dence of the critical temperature on the concentration of ordinary point-like impurities is found. 
Its nonuniversality could serve as a qualitative measure of the two-band character of ferromagnetic 
superconductors. The problem of the upper critical field determination is also discussed. 
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The extension of the Bardeen-Cooper-Schrieffer (BCS) theory of superconductivity to the case with two bands of 
itinerant electrons was developed soon after the appearance of the BCS theory [1]. This theory has been the subject of 
renewed interest following the recent discovery of MgB 2 - the first superconducting material where the existence of two 
energy gaps has been unambiguously demonstrated by thermodynamic and spectroscopic measurements [2] . Certainly, 
there are many other superconducting compounds where multiband effects are less pronounced and experimentally 
invisible because the Cooper pairing occurs mostly in one band of the itinerant electrons or holes. On the other hand, 
there is a whole class of superconductors where two-band (or more generally multiple band) superconductivity is an 
inherent property: the so called ferromagnetic superconductors where the different bands with spin "up" and spin 
"down" electrons are always present. UGe2 [3,4], ZrZn 2 [5], and URhGe [6] are recent examples of such materials 
where the superconducting states are expected to be spin triplet in order to avoid the large depairing influence of the 
■ exchange field due to ferromagnetism. 

The symmetry classification of the superconducting states for itinerant ferromagnetic spin-triplet superconductors 
has been proposed recently by several authors [7-9]. At the same time, a general Gor'kov-type mathematical descrip- 
tion of the multi-band superconductivity in a ferromagnetic metal with a triplet pairing has not been developed yet. 
\Q • The principal goal of this article is to present such a description for a two-band ferromagnetic metal with an anisotropic 
spectrum of quasiparticles and a general form of the pairing interaction. Being useful for a concrete calculation with 
a particular form of the spectrum and the pairing interaction, this approach allows to solve in general terms several 
problems typical for the superconductivity theory, such as: the critical temperature determination, the derivation 
of the Ginzburg-Landau equations, the suppression of the superconductivity by impurities, the upper critical field 
calculation. 

We begin with the general form of the order parameter and the pairing interaction in a two-band itinerant ferro- 
• magnet. Then the Gor'kov equations will be written. That permits to calculate the spectrum of the quasiparticles 
and the critical temperature, and to derive the system of the coupled equations for the order parameters from two 
bands. Then, a law of suppression of the critical temperature by point-like nonmagnetic impurities is found. Its 
characteristic nonuniversal behavior can serve as a qualitative measure of the two-band character of ferromagnetic 
^ , superconductors. The problem of the upper critical field determination is finally discussed. 



II. FERROMAGNETIC SUPERCONDUCTORS WITH TRIPLET PAIRING 

A. Two-band superconductivity 

For a triplet superconductor the order parameter is written as [10] 



A Q/J (R,k)=(£j M = ( d r (R)k)<7 )^ 
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= { -d x (R,k) + id v (R,k) d z (R,k) 
\ d z (R,k) d x (R,k) + id y (R,k) 

where er = (a x ,a y ,a z ) are the Pauli matrices. The superconducting states d r (R, k) with different critical tempera- 
tures in the ferromagnetic crystals are classified in accordance with the irreducible co-representations T of the magnetic 
group M of the crystal [7-9]. All the co-representations in ferromagnets with orthorombic and cubic symmetries are 
one-dimensional. However, they obey multicomponent order parameters determined through the coordinate depen- 
dent pairing amplitudes: one per each band populated by electrons with spins "up" or "down" and one per each pair 
of the bands with the opposite spins (zero spin projection states). Owing to the big difference in the Fermi momenta, 
the pairing of electrons from the different bands is negligibly small. Hence, we shall neglect the pairing amplitude 
with zero spin projection; in another words A = d z (R, k) = will be taken throughout the paper. Also, we limit 
ourself to the consideration of two band ferromagnetic superconductors with a strong spin-orbital coupling, when the 
two-component order parameter has the form 

d r (R, k) = 1 [-(x + iy)Af (R, k) + (x - l y)A l (R, k)] (2) 

which was first pointed out in the paper [11]. Here x, y are the unit vectors of the spin (or, more exactly, of the 
pseudospin [10]) coordinate system pinned to the crystal axes. 

A T (R,k) = -^(R)/_(k), A x (R,k) = »fc(R)/ + (k). (3) 

The functions /±(k) = / x (k) ± if y (k) and the projections fi(k), with i = x,y are odd functions of the momentum 
directions of the pairing particles on the Fermi surface. The general forms of these functions for the different co- 
representations in ferromagnetic superconductors with orthorhombic and cubic symmetries are listed in the paper [7] . 
For instance, in the case of the A\ representation in orthorombic crystal, they are 

f x (k) = k xUl Al + ik y u 2 Al , f y (k) = k y u 3 Al + ik x u A A \ (4) 

where u\ Al , . . . are real functions of k x 2 , k y 2 ,k z 2 . The simple consequence of this is that the nodes in the quasiparticle 
spectrum of the superconducting A-states in orthorombic ferromagnets dictated only by the symmetry are the nodes 
lying on the northern and southern poles of the Fermi surface k x = k y = 0. On the contrary, for the B-states they 
are on the line of the equator k z = 0. 

The order parameter amplitudes rji (R) and r/2 (R) (which are coordinate dependent and complex) have been con- 
sidered (and discussed) in the paper [7] as being equal. This is in general not true. However, even in the general case, 
they are not completely independent: 

»yi(R) - MR)| e ^ R >, ry 2 (R) = MR)^™**. (5) 

Although being different by their modulous, they have the same phase with an accuracy ±7r. The latter property 
guarantees the consistency of the transformation of both parts of the order parameter under the time reversal. 
The BCS Hamiltonian in a two-band ferromagnet with a triplet pairing is 

H= ( k Mk>L a k'a + ^ Yl Vr «/3( k ' k ') a l k +q/2,a a k+q/2,a a k'+q/2,/3 a -k'+q/2,/3, (6) 

k,k',a k,k',q,a,/3 

where the band indices a and (3 are (f , [) or (1, 2), 

h a =i a - ^sg Q H ext /2 + U(r) - e F (7) 

are one particle band energy operators, the functions e a (including the exchange splitting) and the g Q -factor depend 
on the gauge invariant operator — iV + (e/c)A(r) and on the crystallographic directions. In the simplest case of 
isotropic bands without a spin-orbital coupling, gi^ = ±2H oxt /i? xt- U(r) is an impurity potential, H cxt is an 
external magnetic field, 

V x A = B = H oxt + 4ttM, (8) 

and M is the magnetic moment of the ferromagnet. 
The pairing potential interaction is expanded as 

V af3 (k, k') = -V aP v a (k)^*(k'), (9) 
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where 



p T (k) = -/_(k), ^(k) = /+(k). (10) 

It contains four different interaction terms which can be typically classified into two groups corresponding to: (i) a 
pairing between electrons with the same spin polarization (intraband interaction) and (ii) the interband scattering 
terms with V^i = describing the transitions of the pair electron from one sheet of the Fermi surface to the other 
sheet by reversing the pair spin orientation with the help of the spin-orbit coupling. 

When the interband scattering is negligible V^i = — 0, the pairing of the electrons occurs first only in one of 
the sheets of the Fermi surface like in the A\ phase of 3 He. In general, the superconductivity in each band is not 
independent. 

B. Gor'kov equations 

We want to determine the Green's functions of ferromagnetic superconductors in the absence of external pertur- 
bations and impurity scattering. Even under these simple conditions, the system is not spatially uniform due to the 
inherent presence of AttM. If we neglect AitM , the system is spatially uniform. Then, we can write the Gor'kov 
equations in the form 

(iu n -£ kQ )G„(k,w„) + A Q (k)i^(k,w„) = 1 

(iu„ + £ka) i^(k, LO n ) + At(k)G a (k,O = (11) 

where £ka = £k« — £f and u n = irT(2n+ 1) are Matsubara frequencies. The equations for each band are only coupled 
through the order parameter given by the self-consistency condition 

A a (k) = -Tj2J2 E V a , (k,k')F (k',uj n ). (12) 

" k' /3=T,| 



The superconductor Green's functions are 



G «( k '-™)=-Sr# (13) 



F a Q.,u n )= A ; ( ^ , (w) 



where i?k,a = \J £ka + l^a(k)| 2 - Obviously, the superconductivity in ferromagnetic superconductors is non-unitary. 

The Gor'kov equations taking into consideration the magnetic moment 4wM, an external field and non-magnetic 
point-like impurities can be easily written according to the general procedure described in [10]. We shall not overload 
the paper by this and just write the self-consistency equations near the superconducting transition. 

C. The order parameter equations near the superconducting transition 

This system consists of two equations for the order parameter components with spin polarizations "up" and "down" , 

A a (R,r) = -Tj2 J dv'V a ^{v,v')Gp{v',ujJ)G {v',-ujJ) 

n,/3 

x cxp [tr'D(R)] {A/3 (R, r') + E^i&R)} , (15) 
and two equations for the impurity self-energy components 

S Q «,R) J drG Q (r,cD„ tt )G Q (r,-^ l tt )exp^rD(R)] 

x{A Q (R,r) + £ a (w°,R)}, (16) 

where Cj" = u> n + signcj„/2r a , and r a is the quasi-particle mean free time in the different bands. These mean free 
times are related in the Born approximation to the impurity concentration m through 
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— = TrniTVoaU^, (17) 

with u a the amplitude of the impurity scattering, and Nq u the density of electronic states in each band. 
The operator of covariant differentiation is 



The normal metal electron Green functions are 



G a (r,Q n a ) = y^e^^-^ + MBgp.aHe^)- 1 . (18) 

The order parameter components in different bands are determined in accordance with (3): 

A T (R,r) = - m (R)/_(r), A x (R,r) = »»(R)/ + (r). (19) 



D. The critical temperature T c q 

In the absence of an external field let us first find the critical temperature T c0 of a pure ferromagnetic superconductor 
in the (formally) spatially uniform situation of a negligible ferromagnetic moment M = 0. In this case the anomalous 
impurity self-energy part E a (£5™, R) = 0, and from (15) we obtain the system of equations 

m = (9iVi + .9i2?72)A(T c o), 

m = {92m + 92m) K T co), (20) 

where g\ = Vjj(|/_(k)| 2 7Voj(k)), the angular brackets mean the averaging over the Fermi surface, iVoj(k) is the angular 
dependent density of electronic states at the Fermi surface of the band |- Correspondingly 512 = Vj j^|/+(k)| 2 Afoj.(k)), 
.921 = nt<l/-(k)| 2 AWk)}, and g 2 = V lL (\f + (k)\ 2 N Qi (k)). The function A(T) is 

A(T) = 2^T]T^=ln^, (21) 

In 7 = 0, 577 ... is the Euler constant, e is an energy cutoff. 
Thus, similar to [1] the critical temperature is given by 

T c0 = (2 7 e/7r)cxp(-l/ 5 ), (22) 

where g is defined by the maximum value for which the determinant of the system is zero (20) 

9 = (.91 + .92)/2 + V (.91 - .92)74 + 312321. (23) 
In particular, at 312, 321 *C gi 1 .92, the critical temperature is determined by 

g = max(ffi, g 2 ). (24) 



E. The critical temperature dependence on impurities concentration 

Triplet superconductivity is suppressed by non-magnetic impurities [12]. Moreover, the law of suppression of the 
superconductivity is described by the universal Abrikosov-Gor'kov (AG) dependence [13] 

-^t = *(l + ^l)-*(l) (25) 



2 Aqt) V 2 

valid for any unconventional superconducting state and applicable in particular to a concrete unconventional super- 
conductor independently of the pressure [10]. Here ^ is the digamma function. The variable t — T c /T c0 is the ratio of 
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the critical temperature of the superconductor with a given concentration of impurities n\ to the critical temperature 
of the clean superconductor, and x — n\jn\ c = t c /t is the ratio of the impurity concentration in the superconductor 
to the critical impurity concentration destroying superconductivity, or the inverse ratio of the corresponding mean 
free particle lifetimes. The critical mean free time is given by r c = j/nT c0 . This dependence has been demonstrated 
(although with some dispersion of the experimental points) for the triplet superconductor Sr 2 Ru04 [14]. 

Deviations from the universality of the AG law can be caused by the anisotropy of the scattering which takes place 
in the presence of extended imperfections in the crystal. Such a modification of the theory applied to UPt3 has been 
considered previously [15]. However, a complete experimental investigation of the suppression of superconductivity 
by impurities in this unconventional superconductor, in particular the study of the universality of the behavior, has 
not been performed. 

The nonuniversality of the suppression of superconductivity can also be caused by any inelastic scattering mechanism 
by impurities with internal degrees of freedom of magnetic or nonmagnetic origin. For the simplest discussion of this, 
see [16]. 

Finally, universality is certainly not expected in multiband superconductors. Theories for this case have been 
developed with regard to the unconventional superconductivity in Sr2Ru04 (p-wave, two-band two-dimensional model 
[17]) and conventional superconductivity in MgB 2 (anisotropic scattering two-band model [18]). 

A simple modification of the universal AG law for the suppression of the superconductivity by impurities in a 
two-band ferromagnetic superconductor is derived here. Our consideration is limited to the simplest case of scattering 
by ordinary point-like impurities. Then, due to spin conservation, one can neglect interband quasi-particle scattering 
and take into account only the intraband quasi-particle scattering on impurities. At finite impurity concentration the 
system of equations similar to (20) is: 

771 = 5iAi(T)?7i + gi 2 A 2 (T)r72, 

m = 92i Ai(T)fji +S2A 2 (T)fj2 ) (26) 

where 

A,, 2 (T) - * (I) - * (I + j-Ly) + 1» ^ + W ,27) 

Hence, the critical temperature is determined from the equation 

( 9l Ai(T) - 1) ( ff2 A 2 (T) - 1) - g 12 g 21 Ai(T)A 2 (T) = 0. (28) 

In particular, at t/i 2 , gn <C gi, g 2 the critical temperature is determined by the max(T c i, T c2 ) of the solutions of 
the equations 



h £- , G + isy- , (5H-* (r - ) - < 3 »» 

Let us accept for determination that g\ > <? 2 . Hence, the maximal critical temperature in the absence of impurities 
is defined by l/gi — A(T c0 ). Then, at small impurity concentrations the solutions of (29) and (30) are the linear 
functions of impurities concentration: 

7T 

Tci = T c0 - - — , (31) 



T c2 = T c0 -- + - - -!L. (32) 
92 gi 8t 2 

These lines can in principle intersect each other. As a result, an upturn in the critical temperature dependence on 
impurity concentration T c {m) appears. Such a type of deviation of the T c (m) dependence from the AG-law presents 
the direct manifestation of the two-band character of the superconductivity. On the other hand, an absence of strong 
deviations from the universal one-band curve if found experimentally in a ferromagnetic superconductor would mean 
that the superconductivity is developed in one-band with only electrons with " up" spins paired and the " down" spin 
electrons leave normal (or vice versa). 
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Another specific feature of the ferromagnetic superconductors is that even in the absence of an external magnetic field 
the exchange field H cx ~ E ox / hb acting on the electron spins in a ferromagnet produces in addition an electromagnetic 
field AirM ~ 47t/zb&f 3 acting via the electronic charges on the orbital motion of electrons, and suppressing the 
superconductivity [19]. Hence, the actual critical temperature in ferromagnetic superconductors is always smaller 
by the value ~ AnM /H c2 (T = 0) relative to the (imaginary) ferromagnetic superconductor without AttM. The 
upper critical field H c2 is also purity dependent. That is why the impurity concentration dependence of the actual 
T c in a ferromagnetic superconductor might be determined not only directly by the suppression of superconducting 
correlations by the impurity scattering as in any nonconventional superconductor but also indirectly through the 
supression of H c2 . In fact the second indirect mechanism has a negligible influence because the ratio AttM / H C 2(T = 0) 
is of the order of 10~ 2 for superconductors with an upper critical field of the order of several Teslas. 

Thus, the problem of the determination of the critical temperature in a superconducting ferromagnet is at bottom 
the problem of the determination of the upper critical field in a single domain ferromagnet. 



F. The upper critical field 

The equations for the determination of the upper critical field at least near T c are easily derived from the system 
(15)- (16). Keeping only the lowest order gradient terms we have 

A Q (R,r) =-Tj2 J rfr'V/ Q , /3 (r,r')G (3 (r',^/)G (3 (r',-(i/) 

n,/3 

x { (l - (r'D(R)) 2 /2) A P (R, r') + (ir'D(R)) E^dtf, R)} , (33) 

and 

R) -n iM 2 J dvG a {r,Cj n a )G a (v,^ n a ) 

x {(irD(R)) A a (R, r) + S a (w n a , R)} . (34) 

Finding E Q (d>", R) from the last equation and substituting it into (33), we obtain after all the necessary integrations 
the pair of the Ginzburg-Landau equations for the two components of the order parameter 

T]i = Vii&iT]i + V n a 2 ri 2 , 

V2 = V n a 1 r] 1 + V il a2m, (35) 
where the operator ct\ consists of the previously determined homogeneous part and of the second order gradient terms 

&i - (|/_(k)| 2 iVo T (k))A 1 (T) - KfijDiDj. (36) 

The gradient terms coefficients are 

—j 1 i 

K W = (|/_(k)| 2 A 0T (k)^ Tj (k)^ b (k)}— J2 T~yf~ 



2 „Vo l- T «l 3 



+ (/_(k)JVoT(k)^ Ti (k))(/-*(k)JVo T (k)^ Tj -(k)) ^ r " i " T2 Yl -4^2- ^ 

/ Ill < m 



2 UJn^n 
n>0 



The operator a 2 is obtained from (36)-(37) by the natural substitutions 1 — > 2, f — > l> + — y — ■ 

Now, the problem of the upper critical field finding is just the resolution of the two coupled equations (35). There 
are a lot of different situations depending on the crystal symmetry, the direction of spontaneous magnetization and 
on the external field orientation. The simplest case is when the external magnetic field is parallel or antiparallel to 
the easy magnetization axis. If this latter coincides with the 4-th order symmetry axis in the cubic crystal as it is the 
case in ZrZn2, then the gradient terms in the perpendicular plane are isotropic and described by the two constants 
Kfij — K^Sij and K^j = K±5ij. This case formally corresponds to the problem of the determination of the upper 
critical field parallel to the c-direction in the two-band hexagonal superconductor MgB 2 solved in [20]. Then, the 
linearized Ginzburg-Landau equations describe a system of two coupled oscillators and have their solution in the form 
rji = cifo(x) and r/ 2 = 02/0(2;), where fo(x) = cxp(— hx 2 /2), and h is related to the upper critical field by means of 
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\H c2 ±4.M\ = ^°, (38) 

where <£>o is the flux quantum. 

Let us for simplicity limit ourself to the impurityless case. Then, r = 1 — T/T c0 and the equation for the determi- 
nation of the upper critical field is 

[9i(r + A(T c0 )) + V u K } h - 1] [g 2 (r + A(T c0 )) + V^h - 1] 
- [.9i 2 (r + A(T c0 )) + V n K,h - I] [g 21 (r + A(T c0 )) + V^h - 1 ] = 0. (39) 

This is a simple square equation and, as before, if we consider the case 312,321 <C gi, gi and gi > g 2 , then we obtain 
the two following roots 

hi(r) = Tr~r7~i (40) 



^ = ^( T + f-f)- (41) 

This two lines can in principle intersect each other. As a result, an upturn in the temperature dependence of the 
upper critical field appears. 

In the most anisotropic situation such as in the orthorombic crystals UGe 2 and URhGe, all the coefficients K^ xx , 
K^yy, K^ xx and Ki yy are different, even for an external field direction parallel or antiparallel to the easy magnetization 
axis. Then, our system of equations can be solved following a variational approach developed in [20]. Again, an upturn 
in the h(r) dependence can be possible. 

The comparison with an experiment masked by the presence of many ferromagnetic domains will be not easy. 
Monodomain measurements are possible in high enough fields. To work in this region one can easily obtain the fourth 
order gradient terms contributing to the Ginzburg-Landau equations. However, the problem of the theoretical deter- 
mination of the upper critical field at arbitrary temperature has the same principal difficulty as in any conventional 
anisotropic superconductor [21]. 



III. CONCLUSION 



Ferromagnetic superconductors are in general multiband metals. The two-band description of ferromagnetic su- 
perconductors with a triplet pairing developed in this paper presents the simplest model applicable to this type of 
material. One-band superconductivity in these superconductors arises only at a negligibly small spin-orbit coupling, 
as it is the case for the A\— phase of 3 He. We studied the dependence of the critical temperature T c on the concentra- 
tion of ordinary point-like impurities in the framework of a two-band weak coupling BCS theory. We demonstrated 
that the non- universal T c (x) dependence could serve as a qualitative measure of the two band character of the super- 
conductivity in ferromagnetic superconductors. Also, the general equations for the determination of the upper critical 
field at arbitrary temperature and impurity concentration were derived. The solution of these equations near the 
critical temperature was found in the simplest case of a cubic crystalline symmetry for the field orientation parallel 
or antiparallel to the 4-th order symmetry axis. 
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